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ABSTRACT 


The sensitivity of the XR-3 roll behavior in calm 
sea on sidewall parameters and force and moment 
calculations is investigated with the Loads and 
Motions compu ter program and compared with 
experimentally aeasur2d data. Propulsion and rudder 
subroutines and added mass computation are reviewed 
and modified. Recommendations for improved simulation 


of the XR-3 roll behavior are given. 
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A. Background 


In the past several years the Surface Effect Ship (SES) 
has received increased attention in the United States Navy 
and detailed investigations of its sea performance have been 
carried out. The two major categories of SES are the Air 
Cushion Vehicle (ACV) and the Captured Air Bubble Ship 


(CAB). 1 \B-type of craft has become of primary interest 
to the States Navy and two of the constructed 
testc e the Bell Aerospace Systems 100-B and the XR-3 


of appLoximately 100 and 3 tons of displacement, 
respectively. In 1976 the U.S. Navy Surface Effect Ship SES 
100-B of the CAB-type established a world speed record for 
surface type ships of 89.4 knots. In April 1976 the SzS 
100-B launched a missile at 60 knots. The launch was 
successful and the missile hit its target five miles away. 


In order to investigate the characteristics of the 
CAB-type SES under any design and operating conditions 
without the costly need of modifying the actual craft, a 
Loads and Motions digital computer program was developed by 
Oceanics, Inc. [Ref. 1]. In 1970 the XR-3 testcraft was 
delivered to the Naval Postgraduate School and Leo and 
Boncal converted the 100-B Loads and Motions program to 
represent the XR-3. Since there were substantial design 
differences between these two ships programming 


modifications were required in certain subroutines. 


B. Objectives 


The purpose of this thesis is to investigate one aspect 
of saf2 maneuverability of these high speed ships, i.e. 
their roll behavior in a turn maneuver. To fulfill this 
objective the sensitivity of the simulated XR-3 roll 
behavior in turn maneuvers as it is represented by the Loads 
and Motions Program is investigated for its dependence on 
changes in sidewall characteristics. The effect of 
modifications in the added mass computation, sidewall 
parameters and force calculations as well as propulsion and 
rudder input parameters are investigated. An evaluation of 
the computed performance is obtained by comparison with 


experimentally measured steady state roll conditions. 
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II. [INITIAL REMARKS 
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A. PROBLEM OF REPEATABILITY 


It has been this author's experience as well as all 
previous investigators at the Naval Postgraduate School 
using the XR-3 Loads and Motions Program that there have 
been no observed problems related to repeating the results 
for a given input condition. However, during the author's 
first studies of the Loads and Motions Program as listed in 
Menzel's thesis [{Ref. 2] an attempt was made to duplicate 
some of the simulation runs given in that reference work. It 
was found necessary as a first step to restore the 
Simulaticn program on the IBM 360/67 computer from a card 
source onto a disk.! Using the input conditions stated in 
Ref. 2 it was impossible to obtain identical time 
histories, e.g. for turns. Bach time subroutine SIDEWALL was 
read into the computer together with the data input deck, 
thereby overriding the stored program, an error message in 
the printed output indicated that a division by zero 
occurred in this subroutine. The error message was not 
generated when the stored program which contained an 
identical version of subroutine SIDEWALL was called. A 
closer look revealed that in the SIDEWALL-version given in 
Ref. 2 the variable PBAR (plenum pressure) was neither 


1 It should also be noted that results obtained in Ref. 2 
used the IBM 360/67 Release 20.6 system while the results 
obtained in this work used the 21.8B Release with HASP 
installed. 
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defined nor transferred by a COMMON-statement and therefore 
a default value of zero was used during the computation of 
PBHEAD (see Appendix C, SDWL 0440). The missing statement 


PBAR = PB - PINF 


Was inserted to subroutine SIDEWALL. With this correction 
the new results were slightly closer to those given in 
Ref. 2. 


The table given below compares output valu2s for roll 
and pitch angle in a turn at 20 kn with constant thrust and 


a 35 degrees rudder step input under calm sea conditions. 


TABLE I 


without with relative 
PBAR-card PBAR-card Ref. 2 deviation 


Roll angle (deg) 
first peak 7.44 Teod Tate 2.6 % 
avrg at 20 sec 2265 pas beh 3.1 11.9 


ay 


Pitch angle (deg) 
first peak 1.34 1.33 a 9.3 
avrg at 20 sec 0.99 0.94 0.8 14.9 


a 


wn 


Although the differences between columns 2 and 3 in TABLE I 
Might be considered +9 be small as far as magnitudes are 
concerned , the relative deviation is up to 15 % as shown in 
column 4. But what is more important is the fact that the 
roll and pitch angle responses for the uncorrected (plots 1 
and 2 in Appendix A) as well as for the corrected program 
(plots 3 and 4) show an unstable craft behavior for t>40 sec 
where the pitch angle increases rapidly while the roll angle 


decreases rapidly approaching the zero degree value. This 


unstable condition is probably due to the 35 degrees rudder 
step input at 20 kn, as it is used in Menzel's study, which 
is physically unreal since it could never be introduced to 
the XR-3 testcraft in an actual run at that speed. 
Therefore rudder angles of up to 15 degrees introduced at a 
rate of 5 deg/sec will be used throughout this work, as it 
is also done on the testcraft. 


Also, regarding plots 1-4 again and not considering 
whether they really reflect the actual craft behavior, it 
should be noted that a short time history ( ¢<20 sec in 
Ref. 2) could possibly result in wrong conclusions because 
the unstable condition is not evident at that instant of 


time. Therefore, time histories of up to 40 or 50 seconds 


will be shown throughout this study. 


B. STEADY STATE CONDITIONS 


Due to the change in subroutine SIDEWALL as mentioned in 
the preceding section the steady state conditions for 
straight runs in calm water have been reestablished for 
various speeds and are listed in Table II. 


Taple If 


Speed Draft Pitch Plenua Thrust 
angle pressure 

(kn) (in) (deg) (psf) (1b) 
10.0 8.17 1.62 23.93 400.61 
U2a5 1.03 1.11 24.84 349 .42 
15.0 6.74 0.84 24.84 335391 
17.5 6.41 0.63 24.84 346.45 
20.0 6.12 0.48 24.84 373.44 
22.5 5.87 0.36 24.84 411.53 
25.0 5.66 0.29 24.84 453.62 
27.5 5.48 0.25 24.84 513.31 
30.0 5.34 0.26 24.84 574 .43 


These steady state conditions have been established by 
first executing the XR-3 Loads and Motions Program with 
constant speed for 40 seconds and then repeating the run 
keeping the evaluated thrust constant. A comparison of 
these steady state values with those given in Ref. 2. shows 
differences especially in the lower and upper speed range of 
again up to 15 % the cause of which could not be suspected 
anywhere else since the simulation program has not been 


changed since last used by Menzel. 


C. SIGN CONVENTIONS 


The sign conventions used in the XR-3 Loads and Motions 
Program are identical to those used in Report 71-84 by 
Oceanics Incorporated (Ref. 1]. A right handed coordinate 
system is applied to the craft with the positive X-axis 
being measured forward and the lateral Y-axis being measured 
positive to starboard. The vertical Z-axis is measured 
positive downward. Identical signs are used for respective 
forces along those axes. 


The signs of the angles are defined in the following 


Manner : 


pitch angle being positive upward 
(boat noses up) 


roll angle being positive to starboard 
(boat rolls to starboard) 


yaw angle being positive to starboard 
(boat turns to starboard) 


rudder angle being positive to port 
(Tight rudder, boat turns to port) . 


zero pitch and roll angle are referenced to the {X-¥ 


plane parallel to the water surface. 


EFFECTS 
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A. ADDED MASS EFFECT 


1. Theory @ 


A basic element in representing the hydrodynamic 
forces by use of slender body theory are the two-dimensional 
sectional values of lateral added mass A22 of the sidewalls, 
which are necessary for the determination of the lateral 
forces, as well as the two-dimensional sectional vertical 
added mass A33 which is also used in determining the roll 
moment. These two-dimensional added mass values are given in 
the ‘surface Effect Ships Aero/Hydrodynamics Technology 
Design Manual' [Ref. 3] as 


A22 C *RHO*PI*D*D/2.0 


rf 


and 
A33 = C *RHO*B*¥3/8.0 
v 


where C is the lateral added mass coefficient 
C is the vertical added nass coefficient 
v 


RHO is the specific density of sea water 
Dis the local draft 


and B is the width at the ical sidewall waterline. 


The = and C values originally selected in Ref. 1 were 
v 
those corresponding to high speed with di and C =1.0. 
v 
But in order to account for variations in sidewall shape, 


and their influence on the effective lateral added mass, the 
results of the research work described in Ref. 3 led to an 
average value of as This new coefficient improved the 


agreement between theoretical and experimental data as 


stated in Ref. 3 . In the present XR-3 Loads and Motions 
Program a value of sree is used in accordance with Ref. 1 


and the effect of letting ct ge will be investigated in 


later parts of this thesis. 


Considering the equation for the vertical added mass 
A33, a major difference has been found between that one 
given in Ref. 3 and the one used in the XR-3 Loads and 
Motions Program (e.g. Ref. 2), where in accordance with 
Ref. 1 (1971) there was an additional PI-factor in the A33s 
equation. The letter 's' indicates that the computation of 
the added mass is done at the craft's stern. Since Ref. 3 
(1976) reflects the experimental and theoretical work being 
done based on Ref. 1, the A33s equations in subroutines 
SIDEWALL and SIDETAB have now been changed to the new form 
given above. 
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2. New Steady State Conditions 
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The change in the simulated XR~3 performance due to 
the reformulated A33s computation can be observed from 
TABLE III giving the new steady state conditions for 
straight runs in calm sea for various speeds which are 
identical for poth (0.4 and 0.8) lateral added mass 


coefficients. 
TABLE III 
Steady state conditions 
{aew_A33s) 
Speed Draft Pitch Plenum Thrust 
angle pressure 
(xn) (in) (deg) (pst) (1b) 
15.0 6.80 0.86 24.84 336.30 
17.5 6.46 0.65 24.84 347.23 
20.0 6.17 0.49 24.84 374.38 
22.5 5.92 0.37 24. 84 412.62 
25.0 5.71 0.30 24. 84 459.90 
Zi §.53 0.27 24. 84 514. 89 
30.0 5.40 0.28 24.84 576.66 


Comparing these steady state values with those previously 
given in TABLE II in Sect. II.8 one finds in general for all 
speeds that 

- draft has increased by about 0.05 in (0.8 %) 


- pitch angle has increased by 0.02 deg (2 to 7 %) 


? - thrust has slightly increased by less than 0.5%. 


3. Effect in a Turn maneuver 


The effect of the change in the formulation of the 
vertical added mass ina Simulated turn maneuver has been 
studied next. After a 5 sec straight run at 20 kn in calm” 
sea a port rudder deflection was introduced at a rate of 
5 deg/sec and then kept fixed at 15 degrees resulting in a 
turn to port. The deadrise forces are computed at the 
transom. In TABLE IV are shown the steady state conditions 
and the roll moments contributed from the bow seal (FKBS), 
stern seal (FKSS), sidewalls (FKSW), rudder (FKRUD), 
propuision system (FKP), aerodynamics (FKAED) and plenum 
pressure (ABPB*PHI*Z) . 


TABLE _ IV 


=== St a ee Se eS Se = 


=== = SS eS = 


old A33s new A33s deviation 

pitch angle (deg) 0.48 0.52 8.3 % 

roll angle (deg) 2.0 1.98 - 1.0 % 

draft (in) 5.86 5.94 1.4 % 

speed (kn) 19.25 19.27 0.1 % 
moments: 

FKBS -391.1 -386.3 - 1.2 % 

FKSS - 3.6 - 3.6 0.0 % 

PKSW 191.0 202.0 5.8 % 

FKRud - 97.5 -110.0 -12.8 % 

FKP 0.2 0.2 0.0 % 

FKAed - 54.4 - 53.5 1.7 % 

ABPB*PHI*Z 355.4 351.9 1.0 % 


= 
+ 


From TABLE IV it can be seen that, due to the 
modified A33s computation, the counteracting roll acments 
from the sidewall (FKSW) and rudder (FKRUD) changed by about 
6 and 13 %, respectively, resulting in only a little effect 
in generating a tendency toward a smaller roll angle which 
is directed out of the turn. Other major contributions to 
the roll angle are from the bow seal (FKBS directed into the 
turn) and from the plenum pressure (ABPB*PHI*Z) directed out 
of the turn. The corresponding piots (plots 5 - 8 for old 
A33s, plots 9 -— 12 for new A33s) show that with the new 
equation for A33s the pitch and roll angle response curves 
are less damped than they were using the former equation for 
vertical added mass calculation as would be expected. 


It should be noted that the roll moment due to 
rudder is fairly large and changes guite significantly while 
the moment contributed by the propulsion is very smail and 
constant. The cause of this will be investigated in 


chapter IV. 


B. DEADRISE ANGLE 


1. Background 


Surface Effect Ships (SES) of the captured air 
bubble type (CAB) have rigid sidewalls immersed into the 
water, therety - together with the bow and stern seal - 
capturing the air in the plenum chamber and reducing air 
leakage. Since the sidewalls of the XR-3 are not uniform in 
cross-sectional shape throughout their length but are curved 
on the outboard side near the bow similar to a single hall 
ship, the expression "sidehulls" would be more appropriate, 
as can be found in modern literature. 


The importance of the correct understanding of the 
effect of the sidewalls on the craft's performance is 
readily seen from a report on SES Research with the XR-1 
testcraft (Ref. 5]. As reported, the first test series were 
carried out successfully with a 45 degrees sidewall deadrise 
angle, which is the angle between the horizontal plane and 
the ship's outer sidewall surface. In 1964, after a 
modification of the craft to 60 degrees deadrise angle and 
articulating seals, an unstable roll condition occurred 
during a turn at 35-37 knots, a maneuver performed many 

= times before. The testcraft heeled out of the turn , nosed 
down due to the retraction of the bow seal, continued its 
outward roll motion and then flipped over. The report, 
however, does not mention the rudder action generating the 
turn. After this accident, the deadrise angle had been 
changed back to 45 degrees. From this experience the 


tf surface effect ship's roll stability should be expected to 


tee 


be sensitive to the deadrise angle and an investigation of 
this sensitivity is reported in the following sections. 


2. Forces due to Deadrise Angle 
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The deadrise angle, as pointed out by Ref. 3 and 4, 
has a major effect on the craft's roll motion and therefore 
is an important design consideration. During a turn the 
craft is desired to heel inboard, thereby minimizing the 
apparent transverse acceleration (a coordinated turn). For 
typical SES designs, it is not possible to achieve perfect 
coordination, but it is possible to achieve nearly flat 
turns if the craft is designed properly. The principal 
factors affecting roll stability at speed are 


* sidewall geometry 


* seal characteristics 


* vertical location of the center of gravity. 


The forces acting at the center of gravity station of the 
starboard sidewall in a port turn are shown in Figure 1. 
The centrifugal force generated acts away from the center of 
the turn. This force must be counteracted by hydrodynamic 
forces which are created by yawing the ship into the turn. 
To an observer aboard the craft, it would appear as if the 
craft is sideslipping. Wave buildup on the outboard sides of 
the sidewalls increases the pressure there, while only small 
pressure changes occur on the inboard side of the opposite 
sidewall generating a small (negligible) vertical force. 
The principal force arises from the outboard side of the 
sidewall. The direction of this force depends on the slope 
of the deadrise angle. If the resultant force passes above 
the center of gravity (solid line in Fig. 1) a restoring 
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moment results, tending to heel the craft into the turn 
thereby improving its roll stability. If the deadrise angle 
is chosen to be larger (e.g. 70 deg, shown dashed), the 
resultant force passes below the center of gravity, 
resulting in a moment tending to heel the craft out of the 


turn. 


Thus, once the SES geometry and the exact vertical 
location of the center of gravity are known, an approximate 
check for stability in a turn can be made. 


a 


Figure 1 - GENERAL ROLL STABILITY CHECK IN A PORT TURN 
(C.F. CENTRIPUGAL FORCE, €.G. CENTER OF GRAVITY) 


3a SPS 


* 


From the construction data entered into the sidewall 


table (subroutine SIDETAB) Figure 2 has been drawn showing 
the cross sections at several selected stations of the 
starboard sidewall. Since in most simulation runs (which 
will be discussed later) draft was Less than 7 inches, the 
deadrise angles for all stations immersing into the water 
have been calculated by straight line approximation and are 
given in TABLE V. The stations are counted from bow 
(station 0) to stern (Station 28) and are not equally spaced 
for station numbers less than 11. The center of gravity has 
been located experimentally at station 14. The average 
deadrise angle of all stations listed is 64.3 degrees. 


TABLE_V 
Deadrise_angles_ (deg) for _all_ stations 
Straight line approximation for 7'' draft 
Station Deadrise angle Station Deadrise angle 

5 67.6 17 62.1 
6 70.8 18 64.0 
r, 57.9 19 65.3 
8 54.4 20 66.7 
9 $3.7 21 68.1 
10 53.7 22 6.96 
11 $3.7 23 70.3 
12 55.6 24 T69 
13 57.3 25 72.4 
14 58.5 26 73.8 
15 59.7 27 76.9 
16 60.9 28 78.8 
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XR-3 STARBOARD SIDEWALL CROSS SECTIONS 


UNITS IN FT.) 


From the previous discussion the surface effect 
ship's stability behavior is expected to be sensitive to tae 
sidewalls' deadrise angle. In the XR-3 Loads and Motions 
Program listed in Ref. 2 the deadrise angle at the transom 
(78.8 deg) had been used for the computation of the deadrise 
forces. This calculation was implemented into subroutine 


SIDEWALL using the following statements: 


FYH (J) =-A22S*U* (V+XSS*R-ZS*P) 

e 

e 

CTNDR=0.0 

IP(DSS.LE.0.0) GO TO 22 

CTNDR=(ES-3B(1)) /DSS 

IF (THETA.LT.0.0) CTNDR=0.39391 
22 CONTINUE 

FZHOLD (J) =FZH (J) 

PZHDRP (J) =PM1*FYH (J) *CTNDR*PROMOI 

PZH (J) =FZH(J) +PZHDRP (J) 

° 


PK= (FZH (2) ~-FZH(1) )*YSW+FKD-FY*ZS 


The projection of the lateral force FYH(J) at each sidewall 
(J=1 for port, J=2 for starboard side) produced the 
respective ccmponent of vertical force FZH, where FZHDRP is 
the projected force and CTINDR is the cotangent of the 
deadrise‘angle. Tke roll angle had been already taken into 
account in the computation of DSS (draft at the stern). The 
factor PRCMO1 simply provided a means to arbitrarily change 
the vertical projected force for studies undertaken in 
Ref. 2. The factor PM1 introduced the necessary sign change 
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in the sidewall force computation being dependent on the 
craft's side (EM1=-1 for port, PM1=+1 for starboard side). 
The roll moment FK was partly determined from the vertical 
forces of both sidewalls. 


Comparing the deadrise angles given in TABLE V it is 
seen that the transom deadrise angle of 78.8 degrees is 
unigue and the largest along each sidewall and no reason 
could be found why just the transom deadrise angle should be 
used for the vertical force computation. This fact rather 
came up when Leo and BoOncal (Ref. 8) scaled down the 
Simulation program of the 100-B testcraft to create the Xz-3 
Loads and Motions Program. For the 100-B testcraft the 
sidewall cross sections are uniformly shaped throughout most 
of its length. This uniformity does not occur for the XR-3, 
as shown in Figure 2. Therefore it has been supposed that 
another deadrise angle which is more representative for all 
angles existing along the XR-3 sidewalls, e.g. the deadrise 
angle at the center of gravity location (station 14), could 
possibly be more appropriate for the vertical force 
computation. To investigate the effect of this supposition, 
turn maneuvers in calm sea at 20 knots and various rudder 
angles have been simulated for both deadrise angles and 
lateral added mass coefficients. 


TABLE VI 


Steady State Roll Angle (deq)_at_20_kn 


Rudder Deadrise force computation at Experimental 
angle transom center of gravity testcraft 
Cc =0.4 Cc =0.4 Cc =0.8 data (Ref.6) 
h h h 
5 0.52 0.32 0.29 0.09 
9 1.07 0.78 0.72 0.28 
12 1.50 1.15 1.11 0.58 


1.98 1.53 unstable 1.36 


Table VI shows the steady state values for the roll 


angle which are compared with those measured experimentally 
in 1974 and documented in Ref. 6. Since the XR-3 Loads and 
Motions Program under investigation represents th2 XR-3 
craft configuration as it existed in 1974, Ref. 6 may serve 
to check whether the results produced by certain changes in 
the simulation program are favorable or not. It is not 
anticipated to exactly match the results obtained from the 
simulated runs to the measured values since for the 
testcraft data the following precisions for the measuring 
devices are stated in Ref. 6 : pitch and roll angle +0.59, 
rudder angle +1.09 and speed +1 kn. From TABLE VI it is 
obvious that the steady state values for roll angle using 


€ £848 as the lateral added mass coefficient and the center 


of gravity deadrise angle for the vertical force computation 
are closer by 25 to 40 % to the measured angles than are the 
corresponding roll angles considering the transom deadrise 
angle. 


Using the center of gravity deadrise angle and the 
lateral added mass coefficient oe as suggested in 
Ref. 3, the agreement in steady state roll angles could be 


improved by another 3 to 10% for rudder angles up to 
12 degrees. This improvement has also heen reported in 
Ref. 3. Comparing the roll angle responses for turns 
generated by a 12 degrees rudder angle to port with ha 


(plots 21 and 22) and ee (plots 23 and 24) it is found 
that the smaller lateral added mass coefficient generates a 


little higher peak roll angle (1.479) than does the larger 
coefficient (1.279). Both responses show the same number of 
oscillations (eleven cycles) until they die out. The larger 


lateral added mass coefficient generates a small negative 


15 


roll angle shortly after the port rudder motion has been 
introduced. However, for a 15 degrees port rudder angle the 
simulaticn program showed an unstable response with 
increasing amplitudes of oscillation for pitch and roll 
angle (see plots 17 and 18). From this the lateral added 
mass coefficient being 0.8 does not seem to be appropriate 
for the XR-3 simulation if both sidewalls are considered for 
the vertical deadrise force computation. The effect of 
using both lateral add2d amass coefficients will be 
considered again in Sections III.D and £. 


From TABLE VII it can be seen that the roli moments 
contributed by the bow seal (FKBS), plenum pressure 
(ABPB*PHI*Z), rudder (FKRUD) and sidewalls (FKSW) chang? the 
most. Comparing the effect due to the change from transom to 
center of gravity deadrise angle it is found that 

- bow seal inward effect decreases 

- plenum pressure outward effect decreases 


- rudder inward effect decreases 


- sidewall outward effect increases 


" 


outward rall angle decreases. 


The time histories for simulated turn maneuvers with 
15 degrees port rudder angle using center of gravity and 
transom deadrise angle (IABLE VII) are shown in Appendix A 
as plots 13-16 and 9-12, respectively. Comparing these 
plots the effect of the change from transom to center of 


gravity deadrise angle and Oh can be observed as 


- pitch angle response being more damped 


- roll angle response being more damped 


- pitch and roll rate responses being more damped 
and having smaller peak values. 


TABLE VII 
Steady State Values for selected Deadrise Angles 
20_kn_ Turn, calm Sea,_15_ deg port rudder ‘angle 


Deadrise Angle 


(degrees) 78.8 58.5 31.6 28.7 
Pitch angle (deq) 0.52 0.45 0.52 0 .64** 
Roll angle (deg) 1.98 1653 0.31* - 0.33* 
Draft (in) 5.94 5.381 6.07 6.52 
Speed (kn) 19.27 19.37 19.1 ailshengs 
FYSW -824.0 =821e9 -810.0 -801.0 
FYRUD 40.0 Sait 39.5 55.8 
FYAED = 19:5 - 20.0 = en = 17.6 
R*V*AM 803.5 809.2 789.6 762.8 
FKES - 386.3 336.4 - 68.4 4166 
FKSS = 66 ye ET = Ono 065 
PKSW 202.0 212.1 175.0 184.0 
FKRUD -110.0 - 89.9 -108.4 =(53.0 
FKAED = 5363 - 56.0 - 52.4 - 46.2 
FKP 0.19 0.14 0.03 - 0.03 
ABPB*PHI *Z 351.9 21385 54.6 Se we 
PYH (P/S) -39/-147 -45/-126 -78/-95 -107/-88 
FYD -637.8 -651.4 -637.3 -607.0 


arbitrarily chosen 


Note :; * = still decreasing 


** = still increasing, not guite steady state 


The reason for the sidewall roll moment and thereby 
the net roll effect not changing more pronounced (TABLE VIT) 
is that not cnly the lateral force FYH must be considered in 
order to determine the force F effecting roll stability 
(Sect. III.B.2, Fig. 1) but also the lateral drag force FYD 
has to be taken into account. Considering the XR-3 geometry 
using the center of gravity deadrise angle, as it is shown 
in solid lines in Pig. 3, the lateral force FYH can be seen 
to generate the vertical projected force FZHDRP and the 
force F, which is directed well above the center of gravity. 
But since FYD has to be added vectorially to the force F and 
its strength being several times larger than FYH (TABLE VII) 
the new resultant force F' determining the craft's roll 
behavior is directed well below the center of gravity 
(although F being directed above it). Thereby an outward 
roll angle is generated. 


Pigure 3 - DEADRISE AND DRAG FORCES IN A PORT TURN 
(CG= CENTER OF GRAVITY, WL= WATER LINE) 


erty me 


Also from Fig. 3 the effect of considering the 
center of gravity versus the transom deadrise angle on the 
force generating the inward or outward roll moment can be 
seen. The center of gravity deadrise angle and its 
corresponding forces are shown in solid or dashed lines with 
PF being directed above the center of gravity. The transom 
deadrise angle, shown in dotted lines, generates a force P 
being directed below the center of gravity producing an 
outward roll moment which gets reinforced by adding 
vectorially FYD. From this the general stability check 
yields a smaller outward roll moment and angle using the 
center of gravity versus the transom deadrise angle. This is 


in agreeazent with the results given in TABLE VI. 


Next, some simulation runs with arbitrary deadrise 
angles chosen such that small positive and negative steady 
state roll angles result are presented in order t9, first, 
find out which deadrise angle would generate the moments 
necessary to Simulate flat turns, and, second, verify the 
trend of the change in the respective moments which has been 
observed for the change from transom to center of gravity 
deadrise angle in TABLE VII. The value of the lateral added 
mass coefficient in these runs was 0.4 as it was in all runs 
listed in TABLE VII. Considering the roll moments over the 
range of deadrise angles from 78.8 to 28.7 degrees, the 
moments due +o the bow seal, stern seal, propulsion and 
plenum pressure change consistantly in magnitude while those 
due to sidewalls, rudder and aerodynamics reach some extreme 
value and then increase again. Regarding the steady state 
roll angles listed in fTASLE VII, it can be seen that the 
arbitrary chosen deadrise angles happen to be nearly 
symmetric about the deadrise angle that would generate zero 
roll angle (flat turn) under these simulation conditions. 
Keeping this symmetry in mind one finds that the values of 
the roll moments due to bow seal (FKBS), stern seal (FKSS), 
propulsion (FKP) and plenum pressure (ABPB*PHI*Z) are also 


symmetric with respect to zero moments. For propulsion and 
plenum pressure the sign of the moments are identical to 
that of the roll angle while for bow and stern seal mements 
the signs are reverse. Prom the results given in TABLE VII 
with positive steady state roll angles Figure 4 has been 
drawn where the dotted line represents the craft's initial 
State, the dashed line represents the craft's steady state 
considering the transom deadrise angle and the solid line 
represents the craft's steady state (smaller roll angle) 
considering the center of gravity deadrise angle for the 
vertical force calculation. Same markings appiy to the 
moment vectors shown from which their change in magnitude 


for both deadrise angles can be seen. 


A 
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Figure 4 - ROLL MOMENTS IN A PORT TURN 
(ROLL MOMENTS AND ANGLE NOT TO SCALE) 


Legend: AED Aerodynamic CG Center of gravity 
ES Bowseal PP Plenum pressure 
SS Sternseal SW Sidewalls 
RUD Rudder P Propulsion 


Deadrise angle: - - - Transom 


Center of gravity 


Since the simulation run with the center of gravity 


deadrise angle being used for the vertical force computation 
gave results that were more favorable (better damped time 
histories and steady state values closer to those measured 
experimentally) this change in vertical force computation 
was implemented into the XR-3 Loads and Motion Program. 
Starting with the deadrise angle given in Table ¥V for draft 
up to seven inches, again straight line approximations for 
the effective deadrise angle at the center of gravity 
location (station 14) for larger drafts have been performed 
resulting in an almost linear relationship. This part of 
subroutine SIDEWALL is now valid for any draft and contains 
the following statements which have been used in this fora 


for all simulation runs to follow in this study: 


DDRANG=0.0 

IF (DS.GT.0.5833) DDRANG=(DS-0.5833) *0.0629 
DRANG=1.021+D DRANG-PM1¥*PHI 

CTNDR=COTAN (DRANG) 

PZHOLD (J) =F2H (J) 

PZHDRP (J) =PM1*FYH (J) *CTNDR*PROMO1 

PZH (J) =FZH(J) +?ZHDRP (J) 

° 


FK= (FZH (2) -FZH(1) )*Y¥SW+FKD-PY¥*ZS 


It should be noted that vertical projected forces 
from both sidewalls are used in computing the roll moments 
in the XR-3 Loads and Motions Program. Section III.E of this 
thesis will be concerned with the effect of changing the 
Simulaticn program to consider the deadrise projected force 
Only from the outward sidewall in a turn maneuver as it is 


discussed in Sect. III.8.2. 


C. CROSS-FLOW DRAG 


In the preceding section the lateral force FYD which is 
due to cross-flow drag has been found to be rather large 
compared to FYH, thereby changing the direction of the 
resultant force significantly. PYD is computed in 
subroutine SIDEWALL by summing the contributions DF of all 


stations: 


DF(I,J) = -~HRHO*CDSW¥*VREL*A3S (VREL) *DSWAV (I,J) 


where DF (I,J) is the lateral drag contributed by 
the j-th element of each (i-th) sidewall 
HRHO is RHO/2. 
VREL is relative velocity 
DELX is the incremental distance used for the 
wetted draft calculation along each sidewall 
DSWAV is the average wetted draft over one 
incremental distance 
CDSW is the cross-flow drag coefficient for one 
Sidewall . 


The cross-flow drag coefficient being used in the present 
XR-3 Loads and Motions Program is the value corresponding to 
that for flcw normal to a long flat plate being CDSW=1.28 . 
Ref. 3 points out that an investigation of cross-flow drag 
forces and coefficients (before the vertical added mass A33 
had been reformulated as discussed in Sect. III.A.1) showed 
better agreement between theoretical and experimental 
results if CDSW=2.90, representing the normal draq 
coefficient for a sharp flat plate, had been used instead. 
After the reformulation of A33 the agreement between 
modified theory and experiment had very wnuch improved, so 


that the former drag coefficient could be used again. Since 


in this thesis work the change in A33s did not effect the 
roll behavior significantly (TABLE IV), some turn maneuvers 
at 20 kn and 15 degrees rudder angie have been simulated for 
various drag coefficients in order to investigate which CDSW 
would better match simulated and experimental steady state 
roll angle. Prom TABLE VII a decrease in roll angle and 
roll moment is desireable and therefore the cross-flow drag 
force FYD should be decreased via reducing CDSN. The steady 
state values are shown in TABLE VIII with the lateral added 


mass coefficient being unchanged eee 


Steady state conditions for various CDS@ 
20_kn_ turn, 15 deg port rudder angle 
CDSW 1629 1.16 ed) Or 7 
Pitch angle 0.43 0.43 0.44 0.46 
Roll angle 1.55 1.48 1.38 1.14 
Draft 5.78 5.80 5.83 5.91 
Speed 19.39 19.36 19.31 19.19 
FKBS -343.0 -330.5 -308.1 -252.8 
PKSS - 2.7 - 2.6 - 3.8 - 2.0 
PKSW 207.8 175.0 127.8 40.0 
PKRUD - 82.2 - 47.9 - 2.3 80.4 
PKAED - 56.6 - 58.7 - 61.6 - 67.6 
PKP 0.14 0.14 0.13 0.11 
ABPB*¥PHI*Z 276.4 264.6 246.5 201.9 
PYH (P,S) =43/-125 -457-125 -48/-124 -57/-121 
FYD -654.9 -616.9 -561.0 -436.0 


Comparison with Ref. 6 (Table VI) shows close agreement 
in roll angle (1.369) for the run with CDSW=1.0 . But these 
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trial runs have been executed with arbitrary smaller drag 
coefficients because an increase in drag coefficient (e.g. 
CDSW=2.0 for sharp flat plate as in Ref. 3) would have 
enlarged the disagreement between simulated and experimental 
data. Also, from plot 29 for CDSW=1.28 and plot 31 for 
CDSW=1.16 it is seen that a smaller drag coefficient results 
in more damping in the roll response, in this case six 
versus nine cycies of distinguishable oscillations. Since 
the choice of CDSW=1.0 cannot be declared to be appropriate 
for the actual shape of the XR-3 sidewalls it will not be 
considered any further in this study. CDSW=1.16 generated a 
4 % change (about 15 % is desired) toward the experimentally 
measured roil angle and its effect will be investigated 
together with the present cross-flow drag coefficient 
CDSW=1.22 in the studies to follow. 


Shown below in figure 5 is the dependence of steady 
state roll angle in turns generated by a 15 degrees 


L 
angle t0 port at 20 kn on the cross-flow drag coefficients 


used in the simulation runs (TABLE VIII). 


Roll angle 


0.7 0.8 0.9 1.0 1.1 1.2 1.3 


Figure 5 - 5S.S. ROLL ANGLE VERSUS CDSW 


D. THRUST MAPPING 


In turn maneuvers as they have been discussed in the 
preceding sections another important point to be considered 
is the change in thrust on both engines. Since the engines 
mounted on the XR-3 provide constant power for a_ given 
throttle setting, the actual thrust delivered during a turn 
Maneuver may vary even though the throttle remains fixed. 
In ail simulation runs shown in this thesis until now 
constant thrust has been used. In Pef. 6 experimentally 
measured reduction in thrust is given for turn maneuvers at 
20 kn and various rudder angles. From this and Ref. 7 
providing more detailed information the thrust-rudder 


dependence was as listed below. 


TABLE Ix 


—<—= SS Se Se Se a ee ee ee ee eS Se I = 


Se SS ee SS ee Se 


Rudder angle total reduction in thrust (Ref. 6) 
(deg) (%) 
- 2.8 
= 3.3 
12 - 4.3 


TABLE IX kas been used for thrust mapping (Block 16 of 
the data input) in the Simulation runs. Due to 
unavailability of more precise data it has been assumed that 
the total reduction in thrust is contributed from both 
engines in equal amounts. But this may not be quite true, 
Since due to the craft's roll angle a difference in 
immersion depth for both engines could result in unequal 


changes in thrust. TABLE X shows the steady state values of 


these simulations considering two cross-flow drag 


coefficients, deadrise force contributions from both 
sidewalls computed at the center of gravity as well as both 


lateral added mass coefficients Shieh? 


TABLE _X 


=e eh ee ee ee ee ee a Se Sh SS Se 


me a ae ee a ee eS es Se = 


= 


Rudder angle Roll angle (deq) Speed (kn) 
(deg) Simulation kef. 6 Simulation Ref. 6 
CDSW=1.28 
0.30(0. 25) 0.09 19.3(19. 2) 19.3 
0.71(0.61) 0.28 18.9(19.0) 19.5 
12 0.99(0.89) 0.58 18.5(18. 6) 19.2 
CDSW=1. 16 
0.28(0.23) 0.09 19.2(19. 2) 19.3 
Ge67(0.57) 0.28 18.9(19.0) 19.5 
12 0.94(0.34) 0.58 18.4(18.6) 19.2 


Comparing the steady state roll angles obtained if 
thrust mapping is used with those listed in TABLE VI (no 
thrust mapping, deadrise force computation at the center of 


gravity), a change in roll angle by 6 to 14 &% for Sart and 


by 14 to 20 & for aici toward the measured angle 


observed. 


Simulation runs still differ from the roll angles 


Ref. 6 


expected to be a little closer to 


mapping, 
is used. 


steady 


By this also the 


State 


steady 


roll 


factor of about two 


these 


state 


angles found 


can be 


in the 


given in 


three but can be 


if proper 


thrust 


i.e. different maps for port and starboard engine, 


velocities which 


‘ 
i 
: 


differ by about 2.5 % can be expected to come in fetter 


agreement, i.e. to increase for smaller roll angles. 


Regarding the roll angle time responses for simulation 
runs using both lateral added mass coefficients (0.4 and 
0.8) and 2 degrees rudder angle the following 
characteristics have been found: the runs using thrust 
Mapping (Table X, plot 33(35)) snow almost the same 
transient behavior as do the runs without thrust mapping 
(Table VI, plot 21(23)) with about ten cycles of oscillation 
before the transients die out. But the runs using thrust 
Mapping need quite a long time (t>45 sec) before they reach 


the steady state roll angle. 


The effect of using a larger lateral added mass 
coefficient (0.8) can be observed in plot 35. Comparing 
this with plot 33, the roil angle response in plot 35 shows 
a Smaller peak roll angle (1.259 versus 1.479) and the 
transients have smaller amplitudes. But both die out after 
ten cycles. This effect has already been noted in Section 
ELZ.B. 8. 


Investigating the effect of changing the coefficient 
CDSW=1.28 to 1.16 by comparison of the roll responses shown 
in plots 33 and 37 or 35 and 39, a significant damping 
effect as has been found previously in Section III.C can not 


be observed here, 
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E. DEADRISE FORCE OF OUTWARD SIDEWALL 


This part of the investigation on how to effect the roll 
behavior or the testcraft simulation is again concerned with 


the deadrise angle. 


In the present XP-3 Loads and Motions Program deadrise 
forces for both sidewalls are computed with the deadrise 
angle appropriately corrected by the roll angle for port and 
starboard sidewalls. The idea behind this, e.g. in a port 
turn, is to consider a pressure buildup (large PYH(S)) on 
the outboard side of the starboard sidewall and a reduction 
in pressure on the outboard side of the port sidewall (small 
PYH(P)). Thereby an upward vertical force FZHDRP (with 
negative sign) is obtained for starboard side, while a 
downward vertical force (positive sign) is computed ‘for port 
side. Both vertical forces are then added to the buoyancy 
force of their side. The forces obtained by this approach 


are shown in Figure 6. 


In Ref. 4, however, a dirferent philosophy in regarding 
the acting forces is presented. as discussed ia 
Section III.B.2 an approximate check on the craft's roll 
stability can be made if the SES geometry and the vertical 
location of its center of gravity are known (Fig. 1). Ina 
port turn the principal force effecting the roll behavior 
arises from the outboard side of the starboard sidewall 
(relative high pressure on the structure there due to wave 
buildup) while there is only a small force (due to small 
pressure change) on the inboard side of the port sidewall 
contributing a roll moment only if there exists a roll 
angle. Now both vertical forces are directed upward. This 


approach and the resulting forces are shown in Figure 7. 


FZHDRP(P) FYH(S) 


Figure 6 - ACTING DEADRISE FORCES (TWO SIDEWALLS) 


FYH(P) 


FYH(S) 


Pigure 7 - ACTING DEADRISE FORCES (ONE SIDEWALL) 


Note : Since the vertical force contributed by the inboard 
Sidewall is rather small, as shown above, it can be 


neglected in the roll stability check (Sect- III.B.2). 


To investigate the effect of this interpretation of the 
acting forces two statements controliing the computation of 
deadrise forces via the sign of the rudder angle RUDANG 
(port rudder is positive) and the sidewali under 
investigation (PM1=-1 for port side) have been added to the 


computation procedure for the vertical projected force : 


DDRANG=0.0 

IF (DS.GT.0.5833) DDRANG=(DS-0.5833) *0.0629 
DRANG=1.021+D DRANG-PM1* PHI 

CTINDR=COTAN (DRANG) 

RUD SIG=SIGN (1.,RUDANG) 

IF (RUDSIG.NE.PM1) CTNDR=PM1*TAN (PHI) 
PZHOLD (J) =FZH (J) 

PZHDRP (J) =PM1*FYH (J) *CTNDR*PROMO1 

PZH (J) =FZH (J) +PZEDRP (J) 


The steady state values for simulated turn maneuvers using 
thrust mapping and center of gravity deadrise angle are 
shown below considering two cross-flow drag coefficients and 


both lateral added mass coefficients wee 


TABLE XI 


<<< SS 


ee a as a a SS eS ee Se = 


Rudder angle Roll angle (deg) Speed (kn) 
(deg) Simulation Ref. 6 Simulation Ref. 6 
CDSW=1.28 
3 0.43(0.52) 0.09 19.4(19.5) 19.8 
9 0.90(1.00) 0.28 19.1(19.4) 1955 
12 1.22(1. 34) 0.61 18.7(19.0) 19.2 
CDSW=1.16 
5 0.41(0.51) 0.09 19.3(19.5) 19.8 
9 0.87(0.98) 0.28 19.1(19. 3) 19.5 
12 1.17(1.30) 0.61 18.6(19.0) 19.2 
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Comparing the steady state values for simulated runs in 
tables xX and XI , the effect of using only the outward 
Sidewall for the deadrise force computation can b2 obdserved 
to give less favorable (little larger) steady state roll 
angles but mcre favorable (larger) steady state velocities. 
This is exactly what has been expected since, as shown in 
Figures 6 and 7, if deadrise force contributions from both 
sidewalls are used, the deadrise force from the port 
Sidewall will reinforce the one frou the starboari sidewall 
toward a less outward roll moment (Fig. 6) while applying 
the philosophy depicted in Fig. 7 the vertical force from 


the port sidewall counteracts the starboard deadrise force. 


Regarding the roli angle responses shown in Appendix A 
for 12 degrees rudder angle with CDSW=1.28 (plots 41 and 43) 
and CDSW=1.16 (plots 45 and 47) and comparing these with the 
corresponding ones from the previous section (plot 33, 35 
and 37, 39) the change in the vertical force computation 
considering the outward sidewall only is seen to result in 
akrout 8 versus 10 cycles of transient oscillations with less 


amplitude. Using 5 las (plots 43 and 47) has a quite 
significant damping effect and makes the roll angle response 


more resemble a step, which is in accordance with the 
experience gained by Ref. 7. Also, comparing these plots 
with plots 35 and 39, the negative roll angle occurring at 
the moment when the rudder has been introduced does not show 
up any more. 
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F. VERTICAL LOCATION OF CENTER OF GRAVITY 


In Section III.B.2 it was discuss2d how to check for the 
testcraft's roll bahavior in a turn provided the SES 
geometry and the location of the center of gravity are 
known. This section will be concerned with the effact of 
different vertical locations of the canter of gravity on the 


Simulated XR-3 roll behavior 


From experimental neasurenent the location of the center 
of gravity of the XR-3 has been determined to be 10.05 ft 
forward of the transom and Leo and Boncal established the 
vertical location at 2.54 ft above the keel on the 
longitudinal center line of the craft. As modifications to 
the testcraft were introduced (e.g. engines and seals) the 
horizontal location of the center of gravity was again 
established experimentally. No such changes were made on the 
vertical location. In order to investigate the sensitivity 
of the steady state roll angle to vertical locations of the 
center of gravity, reductions in height (in accordance with 
the stability check geometry depicted in Figure 1) in 
increments of 0.1 ft were entered into the Loads and Motion 


Pregran. 


Simulation runs with constant thrust and initial speed 
of 20 kn introducing a 15 degrees port rudder angls, 
deadrise force contributions from both sidewalls computed at 
the center of | gravity, the sidewall drag coefficient 
CDSW=1.28 and the lateral added mass coefficient being 0.4 
were executed and the results are shown in Taple XII. From 
the listed data it follows that with a lower vertical 
location of the center of gravity a tendency toward an 


inward roll angle (here actually a smalier outward roll 


angle) can be achieved. As shown, a decrease in vertical 
distance by 0.1 ft resulted in about 0.1 degree decrease for 
both peak and steady state roll angle (almost linearly 
related). From the corresponding plots as indicated in 
Table XII a roll damping effect due to the relocation of the 
center of gravity can hardly be recognized. The number of 
cycles of transient oscillations and the steady state speed 


were about the same for all runs listed in Table XII. 


Roll conditions in 20 kn turn with constant thrust 


=—= = SS = = ee = = 


vertical location of CG experimental 
2.54 2.44 2.34 2.24 testcraft data 
steady state 
roll angle te5S Twhd F232 V222 1.30 
peak 
roll angle 1667 1.53 46 1.435 
peak 
roll rate ZotS 2.07 h98 1.2.9 
cycles of 
transients 8 ) 8 8 
refer to 
plots 13-14 49-50 51-52 53-54 
steady state 
speed 19.4% 1953 19.3 19.3 18.7 


Table XIII shows the effect of changing the vertical 
location of the center of gravity for simulation runs at 
20 kn with 12 degrees port rudder angle and thrust mapping 
applied, deadrise force computation for both sidewalls at 


the center of gravity and the coefficients being 0.4 for the 
lateral added mass and CDSW=1.28 for the sidewall drag. 
These conditions were similar to those used in Table xX. 


* From Table XIII the effect of reducing the vertical location 
of the center of gravity in 0.1 ft increments over a 
reasonable range can be observed to also result in linear 
reductions in maximum and steady state roll angle for the | 
case if thrust mapping is used. Nonlinear reductions appear 
for the peak roll rate. The number of cycles during the | 
transient period (ten) is again not effected py the vertical 
relocation of the center of gravity. Regarding the 
corresponding plots a very slight damping effect can be 
observed for reduced vertical locations. The results given 


in Table XIII are graphically represented in Figure 8 . 


Roll_conditions in 20 kn turn with thrust mapping 
} and_12 degrees port rudder angle 
vertical location of CG experimental 
2-54 2.44 2.34 2.24 testcraft data 
steady state 
roll angle 0.99 0.91 0.894 0.77 0.58 
peak 
roll angle 1.87 1.36 1.25 1.16 
peak 
roll rate 2-48 14492 1.78 1.58 
cycles of 
& transients 10 10 10 10 
refer to 
plots 33-34 55-56 57-58 59-60 


Steady state 
speed 18.5 18.4 18.4 18.4 19.2 


2.24 2.34 2.44 2.54 


a 2 Figure 8 - ROLL CONDITIONS VERSUS 
VERTICAL LOCATION OF CENTER OF GRAVITY (TABLE XIII) 


G. ROLL CAMFING DUE TO VERTICAL WAVE GENERATION 


One part of subroutine SIDEWALL is concerned with the 
roll danping effect which is due to vertical wave generation 
during a roll motion. The computation of the roll damping 
term was developed by Oceanics Inc. and added to the 
program by Menzel [Ref. 2]. The principles used in tne 
development cf this addition were not provided by Oceanics 
Inc., however, the apparent idea is to reduce the computed 
roll moment FK by some value which is dependent on the 
craft's roll rate and its draft, thereby on the vertical 
added mass A33. During the study of the effect of these 
additional terms a discrepancy was noted in the dimensions 


used in the calculations. 


The roll moment FK is calculated by the equation given 
in SDWL1950 and considering the expressions leading to this 
statement it is found that 


dim [PK] = lb. ft 


which is the correct dimension for a moment vector. Now 


considering the dimension of the roll damping term 
dim [PROMO2*YSW*YSW*BC2*P/PT | 
which is subtracted from FK (SDWL2200) it is found 


PRCMO2 and PI are dimensionless 

dim ([¥SW] = ft 

dim [P] = rad/sec 

dim (BC2] = dim (AC2] 
= dim [Fc2 * length] 
= dim [F2 * cos(x) * length] 
= dim [A33 * length] 


dim —RHO * B * B * length] 7 —~-~—-—... 
(Lb. sece /fts) 2. ft. £f .- ft". 


" 


So for the roll damping tern 
dim [YSW * YSW * P * BC2] 
= £t* . (sad/sec) . (lb «. sece/te*): . £t? 
= lb. ft . sec 


which is not the correct dimension for a moment expression. 
Assumi.g that only the given terms may be involved in the 
roll damping calculation, it is possible to come up with the 
proper dimension if the roll rate P were squared but keeping 
its sign, thus getting 


roll damping term = PROMO2*YSW*YSW*BC2*P*ABS (P) /PI . 


Also the vertical added mass should be associated with some 
"velocity squared' - expression which is absent in the 
original version of this part of the program but could be 
generated by squaring the roll rate P. The effect of the 
modified expression will be an increase in the damping 
effect for roll rates P>1.0 which exist only in the initial 
phase of the turn maneuver (maximum roll rate in all 
Simulation runs was about 2.0) . But for most of the run 
length P is less than unity and the supposed modification 
will decrease the damping effect especially when approaching 
steady state (P=0.0) . 


The PROMO2 - term provides a means to arbitrarily set a 

roll damping factor which, as stated in Reference 2, has 

| been found from experiments to be 16.0. Until now 
PROMO 2=1.0 has been used in this thesis work. So the effect 

of changing it to the experimental value is investigated 

next. In the runs to be discussed turn maneuvers at 20 kn 

with a 12 degrees port rudder angle have been simulated with 

deadrise forces being computed at the center of gravity, 


using the experimental damping factor PROMO2=16.0 and the 


coefficients CDSW=1.28 for the sidewall drag and coer’ for 


the lateral added mass. 


1. Thrust is held constant with deadrise force contributions 
from both sidewalls 

Comparing the obtained plots 61-62 with plots 21-22 
(PROMO2=1.0) only a slight damping e2ffect on roll angle 
response can be observed (peak roll angle 1.459 versus 
1.479), but the number of transient oscillations decreases 
from eleven to nine. The peak roll rate decreased from 2.14 
to 2.10 for PROMO2=16.0 and the damping effect is more 
pronounced in the roll rate time history. 


2. Thrust is mapped with deadrise force contributions from 
both sidewalls 

Ccmparing the plots 63-64 with the corresponding ones 

(33-34) for PROMO2=1.0 , again only a slight damping effect 

can be noticed. The peak roll angle is again reduced from 

1.47 to 1.459 and the number of cycies of transient 


oscillations is uneffected. 


If the damping term is changed to include the P*ABS(P) 
expression (plots 65-66) there result eleven versus nine 
cycles of transient oscillations with Little higher 
amplitudes (peak roll angle is again 1.479 as it was for 
PROMO2=1.0) . The roll rat2 time history also shows little 
larger amplitudes (peak roll rate is 2.14 versus 2.10) . 


3. Thrust is mapped with deadrise force contribution from 
the outward Sidewall only and lateral added mass 
coefficient nee 


Comparing these plots (67~68) for PROMO2=16.0 with plots 
43-44 where PROMO2=1.0 again only a slight damping effect 


in roll angle response can be noticed with the number of 


transient oscillations being uneffected. Using the P*ABS(P) 
term the slight damping effect shown in plot 67 is cancelled 
as shown in flot 69. 


The above observations show that the experimentally 
determined roll damping factor PROMO2=16.0 has only a 
negligible damping effect. Considering the roll damping term 
again as it was added to the present XR-3 Loads and Motions 
Program in 1974, there was a 1/PI - factor which actually 
cancelled the PI - factor contained in the former expression 
for the vertical added mass A33S (Section III.A). Since the 
PI - factor already has been eliminated by correcting the 
A33S - expression, there is no further need to cancel it in 
the damping term by the 1/PI - term which therefore also 
could be eliminated. This change is expected to cause an 
increase in sidewall roll damping by a factor of PI and has 


been investigated by an additional run. 


4. Thrust is mapped with deadrise force contribution from 
the outward sidewall only, roil damping term being 
16. O*YSW#YSW*BC2*D 
and lateral added mass coefficient hs 


From plots 71-72 again only a slight damping effect versus 
plots 67-68 can be noticed. 


The investigation described in this section shows that 
the roll damping due to vertical wave generation as included 
in subroutine SIDEWALL is not very effective even if the 
factor PROMO2 in the damping term is increased from 16.0 to 
50.0 . There also is a lack in matching dimensions 
(statement SDWL2200). 


In all simulation runs executed during this study the 
XR-3 Loads and Motions Program as given in Ref. 2 has been 
used except for the modifications mentioned in this work. 
This simulation program contains Forbes! version [Ref. 9] 
of "SUBROUTINE PROP’ which included some revisions of the 
original version given by Leo and Boncal [Ref. 8]. Since in 
TABLE VII (Section I[II.B.4) the roll moment contributed by 
the propulsion system (FKP) is rather small (0.1 to 0.2) an 
additional run using Leo and Boncal's version has been 
executed under the conditions as stated for TABLE VII 
(center of gravity deadrise angle). The steady state values 
are shown below (in parantheses are the values obtained 


using Forbes’ version) : 


pitch angle 0.44(0.45) draft 5.73 ( 5.81) 
roll angle 1.67(1.53) speed 19.39(19.37) 
FYSW -909.4(-821.9) FYRUD <- 546.7( 32.7) 
FYAED - 21.3(- 20.0) R*¥V*¥AM 888.4( 809.2) 
FYP 97.0( 0.0) 

PKBS -348.0(-336.4) FKSS - 3.0(- 2.7) 
PKSW 269.9( 212.1) FKROD 150.3(- 89.9) 
PKAED - 63.5(- 56.0) FKP -305.2( 0.14) 


ABPB*PHI*Z 299.5( 273.5) 


PYD -729.2(-651.4) FYH (P/S) -43/-137(-45/-126) 


Regarding these values drastic differences are found for 
FYP and FYRUD as well as for FKP and FKRUD, though <the 
steady state values for pitch and roll angle, draft and 
speed are close to those obtained using Forbes' version. 
The pitch and roll angle responses are given in plots 25-28. 
Comparing these with the corresponding ones using Forbes! 
version (plots 13-16) it is seen that Leo and Boncal's 
version generates more oscillatory roll and pitch angle 
responses. Therefore a review of subroutine PROP especially 
the computation of the respective forces and moments has 


been carried out and are discussed below. 


The review resulted in force and moment equations which 
are identical to those given by Leo and Boncal {Ref. 8]. 
Ancther point that supports the supposition to use this 
version of subroutine PROP is the fact that it provides a 
reasonable roll moment due to the propulsion system (Forbes! 
version provides almost zero roll moment). In order to test 
the effect of changing subroutine PROP to its original fora, 
Simulation runs 1. and 2. from the preceding section have 
been repeated and their time histories for the roll behavior 
are shown in flots 73-74 and 75-76, respectively. Comparing 
these with the corresponding ones using Forbes' version of 
Subroutine PROP, plots 61-62 and 63-64, it is observed that 
the steady state and peak roll angles are larger by about 
7% and the number of transient oscillations is uneffected 
by this change. 


Next a review of subroutine RUDDER has been carried out 
in order to investigate the cause of the change in sign for 
FYRUD and FKRUD depending on the propulsion subroutine used 
(Forbes' or Leo and Boncal's version) although identical 
RUDDER-subroutines have been used in both cases. The 


important equations to be considered are : 


DSR=Z+ZS-XR*THETA 

ENDFAC=%3. + DSR/(DSR+RSPAN) 

VH=V + XR*R - ZR*P 

EFFANG=RUDANG - ENDFAC*VH/U 
PY=RHO*U*U*RA REA*ENDPAC*RCLB*EFFANG 


FK= -ZR*FY 


The above equations compute 

- draft at the rudder location 

- endfactor depending on the rudder Length below the 
craft's keel 

- sideslip velocity 

- effective rudder angle depending on the craft's forward 
and sideslip velocities 

- lateral force depending on rudder area and effective 
rudder angle 

- roll moment due to rudder . 


The sign of the roll noment FK is dependent on the sign 
of the lateral force FY which depends on the sign of the 
effective rudder angle EFFANG. This in turn depends of 
course on the introduced rudder angle, but the magnitude of 
the second term of this equation (sideslip angle) is 
responsible for the actual sign of EFFANG. VH/U represents 
the tangent of the sideslip angle which ~ applying small 
angle approximation - is subtracted from the introduced 
rudder angle (see Figure 9). Since ENDFAC=1.73 and OJ was 
the same for both runs, a difference in VH which was larger 
in Leo and Boncal's version by about 50 % caused the 
negative sign (versus a positive sign using Forbes' version) 
for the effective rudder angle EFFANG and thereby FY. So 
the only way to introduce a change toward an inward roll 
behavior here (FKRUD was positive) is to first check the 
input value for RSPAN and second to eventually modify the 
craft's engines to have larger immersion into the water. 


With regard to the first item, if RSPAN would increase, 
ENDFAC would decrease as would the magnitudes of EFFANG 
(negative), FY (negative) and FK (positive). 


U 


SSA = VH/U sideslip angle 
ESSA=ENDFAC*VH/U effective sideslip angle 
RUDANG rudder angle 
EFFANG effective rudder angle 
Figure 9 - EFFECTIVE RUDDER ANGLE IN A TURN MANEUVER 


Considering next the fact that the engines of the XkR-3 
' had been replaced in 1976 a recent measurement of the rudder 
and propeller location yielded the following data (original 
input values shown in parantheses) 
XPO = -1.50 (-1.275) 
ZPO -1.333 (-0.604) 
XRO = -1.083 (-1. 125) 
RSPAN= 1.333 ( 1.21 ) 
RAREA= 1.42 ( 0.68 ) 


Two of the runs discussed in the preceding section have 
been repeated using the original roll damping term and the 


new propulsion and rudder input parameters. 


1. Thrust is mapped with deadrise force contributions’ from 
both sidewalls with Roneee A 


Comparing plots 77-78 with the corresponding ones for the 
former rudder and propulsion parameters (plots 63-64) the 
expected effect due to a decrease in FKRUD can be observed 
resulting in smaller peax roll angle (1.349 versus 1.450) 
and a significantly decreased steady state roll angle (0.659 
versus 0.999 at t=50 sec). The transient period of the roll 
angle response in plot 77 shows the same number of transient 
oscillations as plot 63 but the response did not reach 
steady state even after 45 sec . The steady state roll 
angle could be estimated to be about 0.3 to 0.4 degrees. 


2. Same conditions as in the previous run but with deadrise 
force contribution from the outward sidewall only and 


c =0.8 . 
n 


Regarding plots 79-80 versus plots 71-72 a similar effect as 
mentioned pefore is observed. Now the roll angle response 
showS a more oscillatory transient period. The peak roll 
angle reduced from 1.46 to 1.36 degrees while the roll angle 
at t=50 sec changed from 1.34 to 0.97 degrees. For this run 
the steady state roll angle could be estimated to be about 
0.8 degrees since the roll angle response again did not 
reach steady state for the run using the recently determined 
rudder and propulsion parameters. The initial small 
negative roll angle which appeared in the previous run 
(plot 77) at the moment when the rudder motion is 
introduced, does not show up any more. The overall roll 
response shown in plot 79 does not fall off as steeply 
(compared to plot 77) if, in a turn maneuver, deadrise force 


contribution from the outward sidewall only is considered. 


Prom these runs it is obvious that the rudder and 
propulsion parameters o9f the new engines result in a more 
favorable XR-3 roll behavior with a significantly less 
steady state outward roll angle. At the time of this 
sensitivity study Reference 7 pointed out that there 
actually is no roll angle indicated by the measuring devices 
in a port turn at 20 knots and 12 or 15 degrees rudder 
angle. Keeping the device sensitivity of +0.5 degrees in 
mind, it can be concluded that the steady state roll angle 
provided by the present XR-3 Loads and Motions Program 
including the berore mentioned progran nodifications by the 
author, is quite well in agreement witn the actual craft 
data. Further efforts should be undertaken to obtain a 


better damping effect during the transient period. 


Remark ; 

Altnough the new engines extend about 0.5 ft deeper into the 
water than the old engines, the new input value for RSPAN is 
close to the former value. In Reference 1 the denominator of 
the second term in the equation for ENDFAC is defined to be 
the distance from the water surface to the bottom tip of the 
rudder. The equation for ENDFAC used in the present program 
accounts for changes in draft DSR in both the numerator and 
the denominator. The RSPAN-term should then be defined to 
be the distance from the craft's keel to the bottom tip of 
the rudder. The author suspects that formerly there has been 
chosen too large a value for RSPAN in the case of the old 
engines (RSPAN=1.21 ft, probably from water surface to 
bottom tip of the rudder). Since the bottom tip of the 
rudder and the center of the propeller were located at 
almost the same distance from the craft's keel the 
magnitudes of ZPO and RSPAN should be about the same which 
is not true for the former chosen values. From these 


observations the author concludes that until now a too 


Optimistic roll behavior resulted due to a too large value 
for RSPAN. 


The preceding sections of this chapter reflect the 
results of an investigation on how the XR-3 roll behavior in 
turn maneuvers at 20 kn simulated with the Loads and Motions 
Program is affected by certain changes in parameters and 
forces. The results of this sensitivity study are 


summarized as follows : 


* Added mass effect 


A ceformulation or the vertical added mass computation which 
reduced the A33s-value by a factor of 1/PI had only a little 
effect in changing steady state draft, pitch angle and 
reguired thrust for constant speed in straight runs. In 
turn maneuvers this change resulted in less damped responses 
for pitch and roll angle and affected the steady state 
values for pitch angle by +8.3 % , roll angle by -1.0 % and 
draft by +1.4 % . Considering the roll behavior, the effect 
of reducing the A33s-value by a factor of 1/PI in the 
vertical added mass computation was negligible. 


A change in the lateral added mass coefficient ‘eo to 


0.8 had the effect of reducing the amplitudes of the 
transient oscillations and the steady state roll angle by 
0.1 degree toward the experimentally measured value. Since 
in a turn maneuver with 15 degrees rudder angle an unstable 
craft behavior showed up which had not been experienced in 
practice, the value of aa is not realistic if deadrise 


~- 


force contributions from both sidewalls are considered. if 


deadrise force contribution from the outward sidewall only 
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is used in the simulation vorogram the increase in lateral 
added mass coefficient causes an increase in the. steady 
state roll angle by about 0.1 degree. But since now the 
roll angle response is quite significantly damped and 
resembles a step as it does in practice, the choice of the 
larger coefficient (0.8) seems to be appropriate for the 
case of the deadrise force contributed from the outward 


sidewall cnly. 


* Deadrise angle effect 


The deadrise angles over the length of the curved XR-3 
sidewalls are not uniform. Changing the deadrise angle which 
is used in the deadrise force computation from the transom 
(78.89) to the center of gravity (58.59) resulted in more 
damped respenses for pitch and roll motion in turn 
Maneuvers. For all chosen rudder angles generating port 
turns at 20 kn the steady state roll angle changed favorably 
by 23 to 38 % toward the corresponding roll angle measured 
experimentally and documented in Ref. 6. 


* Cross-flow drag coefficient 


Changing the cross-flow drag coefficient (CDSW=1.28 for a 
long flat plate) to lower values, e.g. CDSW=1.00 , the 
testcraft's simulated roll behavior could be forced to match 
the measured steady state values. But since this 
coefficient has been chosen arbitrarily and cannot be 
declared to be appropriate for the actual XR-3 sidewalls 
(CDSW=1.0 means that there is no effective drag coefficient) 
this CDSW-value has not been considered any longer. Another 
Suggested cross-flow drag coefficient CDSW=1.16 [Ref. 3] 
resulted in about 5 % reduction in steady state roli angle 
and a shorter transient period (about 20 &) with smaller 


amplitudes of oscillation than did the runs using CDSW=1.28. 


* Thrust mapping 


Performing thrust mapping in the simulation runs using the 
data provided in Ref. 6 and assuming that both engines’ lose 
egual amounts of thrust in turn maneuvers, th2 simulated 
steady state roll angles could be shown to agree better with 
the measured angles by about another 10 % but the roll angle 
response did not quite reach steady state even after 45 


seconds. 


* Deadrise force from outward sidewall only 


The Simulation runs executed with the deadrise force 
computed for the outward sidewall only and thrust mapping as 
done before, showed less agreement for both ie (0. 3) 


with the corresponding measured steady state roll angles 


than did the runs considering deadrise forces from both 
sidewalls. But for “aes the roll angle response was 


significantly damped and looked more like a step, as was 


experienced in practice. 


* Vertical location of the center of gravity 


The given value for the vertical location of the center of 
gravity (2.54 ft above the keel) might not be quite correct 
due to the difficulties in experimentally determining it. 
The significant effect of the vertical center of gravity 
location on the XR-3 roll behavior is obvious from Figures 1 
and 8 and the discussion in Section III.F. Since it controls 
exceedingly the XR-3 roll behavior it is important that 
future investigators obtain an accurate measurement of the 


C.G. location before taking testcraft data. 


* Subroutine PROP 


A review of the propulsion subroutine revealed that the 
force and moment calculations as given by Leo and Boncal 
(Ref. 8) should be used instead of the version stated by 
Forbes (Ref. 9). This change results ina 7 & larger outward 
steady state roll angle but it provides a reasonable roll 
moment due to the propulsion system. 


* Roll damping due to vertical wave generation 


The roll damping calculation inciuded in subroutine SIDEWALL 
has only a negligible effect on the craft's net roll 
behavior over the range of the damping factor 
1.0<PROMO2<50.0 . Since there also is a discrepancy in the 
dimension of the roll damping term the source of this part 
of the program should be investigated. 


* Suggested areas for further investigation 


During the studies undertaken with the XR-3 Loads and Motion 
Program the author found the following points either in the 
input data from the testcrafr+ or in the program that require 
further investigation : 


s Rudder/thrust angle 
As learned from Ref. 7 the rudder/thrcust motion is directly 
introduced to the port engine while the starboard engine 
follows this via a metal bar connecting both engines. The 
starboard rudder and thrust vector may be off by 2 or 3 
degrees to either side, depending on the direction of the 
rudder/tarust action. Thereby the roll behavior will 
certainly be effected. Since rudder/thrust mapping is 
available in the Loads and Motion Program, it is important 
that accurate data be collected for each engine during a 


test run. 


£ 


e XR-3 weight 

During each test run there is a possibility that the weight 
is not exact since, as learned from Ref. 7, the testcraft 
takes on water in the seals and thereby its weight increases 
causing larger draft which also effects the roll behavior. 
This additional weight can be substantial and attempts 
should be made to measure or estimate the additional weight 
to be used in the Loads and Motion Progran. 


e cCross-flow drag coefficient 
Further efforts should be undertaken to establish the proper 
value of the cross-flow drag coefficient since CDSW=1.28 
might not be quite appropriate for the actual shape of the 
XR-3 sidewalls. 

e Sideslip velocity 
Those parts of the Loads and Motions Program where terms 
contributing to the lateral velocity V are computed should 
be reviewed since a large value of V may result in a change 
in the effective rudder angle from positive to negative. 

In this sensitivity study it has been shown that the 
most significant effect on the simulated XR-3 roll behavior 
are contributed by 

* deadrise angle 

* vertical location of the center of gravity 

* rudder and propeller location . 

Less significant effects are contributed by 


* thrust mapping 


* cross-flow drag coefficient . 
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Por further studies concerned with the XR-3 Loads and 
Motions Program for improved roll behavior representation it 


is recommended to use 


- the reformulated equation for the vertical added mass 
A33 


- deadrise angle at the center of gravity location 


- deadrise force contribution from the outward sidewall 
only 


- subroutine PROP as given by Leo and Boncal and in 
Appendix C of this thesis 


- thrust and rudder mapping. 


f 
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APPENDIX A 


PLOTS 


In the table given below are listed all plots with the 
respective parameters used in the corresponding simulation 


run. The table contains the following abbreviations : 


COEF coefficient CG center of gravity 

CDSW sidewall drag e lateral added mass coef. 
ch4 chapter IV 

DRANG deadrise angle DF16 PROMO2=16.0 

PA pitch angle PR pitch rate 

RA roll angle RR roll rate 

REFTAB refer to table RUDANG rudder angle 

TR Transom 

L-B Leo and Boncal's version of subroutine 'PROP' used 


P*|P| or P*FI P replaced by ... in SDWL2200 (Section III.G) 
ZS~0.1 vertical location of CG is 2.54 - 0.1 (ft) 


PLOT REP RESP RUD DR COEF COEFF THRUST SDWL- REMARK 


NO TAB ANG ANG * CDSW MAP FORCES 

1 1 RA 35 TR 0.4 1.28 no both no 

2 1 PA 35 TR 0.4 1.28 no both ee 
3 1 RA 35 TR 0.4 1.28 no both with 
4 ' 35 TR 0.4 1.28 no both ee 
3 4 RA 15 TR 0.4 1.28 no both 

6 4 FR 15 TR 0.4 1.28 no both old 
7 4 PA 15 TR 0.4 1.28 no both A33s 
8 4 PR 15 TR 0.4 1.28 no both 
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PLCT REF RESP RUD ODOR COEF COEF THRUST SDWL- REMARK 
NO TAB ANG ANG ee CDSW MAP FORCES 
9 4,6 RA 15 TR 0.4 1.28 no both 
10 4,6 RR 15 TR 0.4 1.28 no both 
11 4,6 PA 15 TR 0.4 41.28 no both 
12 4,6 PR 15 TR 0.4 1.28 no both 
13 6 RA 15 CG 0.4 1.28 no both 
14 6 KR 15 CG 0.4 1.28 no both 
15 6 PA 15 CG 0.4 1.28 no both 
16 6 PR 15 CG 0.4 1.28 no both 
17 6 RA 15 CG 0.8 1.28 no both 
18 6 KR 15 CG 0.8 1.28 no both 
19 6 RA 12 TR 0.4 1.28 no both 
20 6 RR 12 TR 0.4 1.28 no both 
21 6 RA 12 CG 0.4 1.28 no both 
22 6 RR 12 CG 0.4 1.28 no both 
23 6 RA 12 CG 0.8 1.28 no both 
24 6 RR 12 CG 0.8 1.28 no both 
25 chy RA 15 CG 0.4 1.28 no both L-B 
26 ch4 RR 15 CG 0.4 1.28 no both L-3 
27 ch4 PA 15 CG 0.4 1.28 no both 
23 ch4 PR 15 CG 0.4 1.28 no both L-B 
29 8 RA 15 CG 0.4 1.28 no both 
30 8 RR 15 CG 0.4 1.28 no both 
31 3 RA 15 CG 0.4 1.16 no both 
32 8 RR 15 CG 0.4 1.16 no both 
: 33 10 RA 12 CG 0.4 1.28 yes both 
34 10 RR 12 CG 0.4 1.28 yes both 
35 10 RA 12 CG 0.8 1.28 yes both 
36 10 RR 12 CG 0.8 1.28 yes both 
37 10 RA 12 CG 0.4 1.16 yes both 
38 10 RR 12 CG 0.4 1.16 yes both 
; 4 39 10 RA 12 CG 0.8 1.16 yes both 
40 10 RR 12 CG 0.8 1.16 yes both 
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42 
43 
44 
45 
46 
47 
48 
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REF 
TAB 


11 
11 
11 
11 
11 
11 
11 
11 
12 
Ue: 
12 
12 
12 
12 
13 
13 
13 
13 
13 
13 
ch4 
ch4 
ch4 
ch4 
ch4 
ch4 
ch4 
chy 
ch4 
ch4 
chy 
ch4 


RESP 


RUD 
ANG 


12 
12 
12 
12 
12 
12 
12 
12 
15 
15 
15 
15 
15 
15 
12 
12 
12 
12 
12 
12 
12 
12 
12 
12 
12 
12 
12 
12 
Zz 
12 
12 
12 


CG 


SE ee ee pee ee neg ge ee ern terre it 
J 7 


COEF 
CDSW 


1.28 
1.28 
1.28 
1.28 
1.16 
1.16 
1.16 
1.16 
1.28 
1.28 
1.28 
1.28 
1.28 
1.28 
1.28 
1.28 
1.28 
1.28 
1.28 
1.28 
1.28 
1.28 
1.28 
1.28 
1.28 
1.28 
1.28 
1.28 
1.23 
1.28 
1.28 
1.28 


THRUST 
MAP 


yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 

no 

no 

no 

no 

no 

no 
yes 
yes 
yes 
yes 
yes 
yes 

no 
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PLOT 63 
for applied parameters see first page of this appendix 
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PLOT 65 
for applied parameters see first page of this appendix 
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PLOT 66 
for applied parameters see first page of this appendix 
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PLOT 70 
for applied parameters see first page of this appendix 
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for applied parameters see first page of 
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PLOT 72 
for applied parameters see first page of this appendix 
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fw for applied parameters see first page of this appendix 
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for applied parameters see first page of this appendix 
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PLOT 76 
for applied parameters see first page of this appendix 
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for applied parameters see first page of this appendix 
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f for applied parameters see first page of this appendix 
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APPENDIX B 


MODIFICATIONS OF THE SIMULATION PROGRAM 


—~ During the course of studies undertaken with the xXR-3 
Loads and Motion Program some statements have been added or 
changed in order to improve the work with the simulation 


program. 
* MAIN program 


The statements 


COMMON /AIR/ PINF,nxHOINF,GAM,IQUIT (MAINO 050) 
IF (MYTIME(DUM) .LT.IQUIT) GO TO 12 (MAINO66 1) 


have been modified and added to the MAIN-program. The 
purpose of the second statement is to start the writing of 
the output values desired if computer time has reached the 
demanded time (jobcard) reduced by IQUIT which is the 
expected amount of time (in 10 * seconds) required to print 
the output. The user may input the desired value of IQUIT as 
the fourth parameter (110-Format, columns 36-45) on the 107 
data card. Default value for IQUIT is 600000 (60 seconds). 


* Subroutine INCON 
The added statements are 

IQUIT=600000 (INCNO 88 1) 
to set the default value and 


READ (5,2041) IQUIT (INCN1455) 
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to read the set value for IQUIT thereby overriding the 
default value. 
* Subroutine RHS 


In Menzel's version [Ref. 2] of the XR-3 Loads and Motions 
Program the entered craft's velocity UO (kn) is changed in 
Subroutine INCON to U (ft/sec) by 


U = U0*1.6889 (INCN4580} 


and in Subroutine RHS transformed back to units in knots 
(RHS 1850) by 


VEL = 0.5925*07 . (RHS 1850) 
From these two eguations follows 
VEL = 0.5925*00*1.6889 


and for JO=20 kn there results VEL = 20.0135 kn . But since 
it is desired that VEL=U0 the statement has been corrected 
to 


VEL = U/1.6889 . 


Statements MAINO840, RHS 1140 have been changed accordingly. 


* Subroutine INTGRL 
To this subrcutine the statement 
CALL COLFIL (INT 1021) 


has been added in order to provide the output values already 
calculated to the user if the minimum time interval allowed 
in the integration process is undergone. Before this 
addition only the warning message 


"DELTA TIME LESS THAN 1.0E-6 - - JOB STOPS! 


appeared and the output values have been lost. 
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